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Abstract 

Correlations of two protons emitted near midrapidity in p + Pb collisions at 450 
GeV/c and S + Pb collisions at 200^1 GeV/c are presented, as measured by the NA44 
Experiment. The correlation effect, which arises as a result of final state interactions 
and Fermi-Dirac statistics, is related to the space-time characteristics of proton emission. 
The measured source sizes are smaller than the size of the target lead nucleus but 
larger than the sizes of the projectiles. A dependence on the collision centrality is 
observed; the source size increases with decreasing impact parameter. Proton source 
sizes near midrapidity appear to be smaller than those of pions in the same interactions. 
Quantitative agreement with the results of RQMD (vl.08) simulations is found for p+Pb 
collisions. For S + Pb collisions the measured correlation effect is somewhat weaker than 
that predicted by the model simulations, implying either a larger source size or larger 
contribution of protons from long-lived particle decays. 

PACS codes: 13.60.Rj; 13.85.-t; 25.75.-q; 25.75.Gz 

Keywords: two-proton correlations; proton-nucleus, nucleus-nucleus; heavy-ion collisions; 

proton source size; CERN SPS 



1 Introduction 

Correlations of identical particles emitted with similar velocities are commonly used as a 
tool to measure the space-time development of the emission process in particle and heavy- 
ion collisions. The technique is most often used for bosons, for which the effect of quantum 
(Bose-Einstein) statistics is the only source of correlations (photons), or plays a dominant 
role (pions and kaons). Other effects, in particular those due to the Coulomb interaction for 
pairs of charged particles, are usually treated as a correction. In the case of protons, final 
state interactions, due to strong and Coulomb forces, typically dominate the effect of quantum 
statistics and determine the form of the correlation function. A negative correlation effect, 
due to Fermi-Dirac statistics and Coulomb repulsion, competes with the positive (attractive) 
correlation due to the strong interaction, giving rise to a characteristic "dip-peak" structure 
(minimum at zero, maximum at ~ 20 MeV/c proton momentum in the pair rest frame) in the 
two-proton correlation function. The height of the correlation peak decreases as the proton 
source size increases. All effects are strongly sensitive to the space-time parameters of the 
emission process [|], ^ 0, |[. 

Correlations of protons emitted with small relative momenta were observed for the first 
time in vr" -|- Xe interactions at 9 GeV/c [0. Currently, two-proton correlations are used 
extensively in low energy (of order 10-1000 MeV/c) heavy-ion physics @, 0, ^, where the 
sizes measured for small proton momenta are comparable to or exceed the size of the larger 
(usually target) nucleus. A similar result was obtained at much higher CERN SPS energies 
in the target fragmentation region, where a two-proton correlation analysis yielded source 
sizes compatible with the sizes of the target nuclei [1^. A decrease of the measured size 
with increasing proton momenta is reported in many papers however, and for very different 
projectile energies [p, ^ ^ |ll|, . Several authors have also reported a dependence on collision 



impact parameter in low energy heavy-ion interactions P, |T3[. The observed trends in the 
measured sizes may be related to differences in the proton emission time |]^, ^, and may also 
reflect the reaction dynamics, which can produce correlations between the momentum and 
position of the emitted particles []T^ . 



In spite of extensive analysis of pion and kaon correlations in relativistic heavy-ion colli- 
sions, data for two-proton correlations are scarce. This paper reports the first measurement 
of proton-proton correlations at midrapidity at CERN SPS energies. The origin of protons 
emitted in relativistic heavy-ion collisions differs from that of lighter particles. Protons are 
constituent parts of the projectile and target nuclei, unlike pions and kaons which are pro- 
duced in the collision. The different dynamics of proton emission may therefore be reflected 
in different source sizes. At the SPS, the expected net baryon-free region at midrapidity for 
the highest collision energies is not yet achieved. There is still significant stopping, with the 
associated rapidity shifts of projectile and target nucleons leading to an increase in the proton 
density at midrapidity |jl5|, |T^. Consequently, the ratio of antiprotons to protons at midra- 
pidity is approximately 0.1 |TB|. The relatively high degree of stopping implies that secondary 



interactions of pions, kaons and other produced particles with nucleons are important, and 
may influence the size of the proton emission region. The relation between the proton source 



size, the dimensions of the colhding nuclei, and the corresponding source sizes for pions and 
kaons at midrapidity for CERN SPS energies is still an open question. In this context, the 
results of the two-proton correlation analysis described here provide complementary informa- 
tion to existing data from two-boson correlations which can be useful in understanding the 
underlying dynamics of the collisions ||T^ |18|, |19], ^, |21|, |2^ . 



The paper presents the results of a proton-proton correlation analysis for collisions of two 
different projectiles, protons and sulphur nuclei, with a lead target. In the first case, the 
projectile consists of a single proton, and the measured pairs contain mainly protons from 
the target or protons produced in the interaction. In the second case, the correlated pairs 
contain both projectile and target participants as well as produced protons. The relative 
role of these different sources will depend on the impact parameter of the collisions, i.e. the 
centrality selection. This dependence is examined by an analysis of the correlation effects as a 
function of charged particle multiplicity in the collision. Finally, the experimental correlation 
functions are compared to predictions obtained from both a Gaussian source model and from 
the RQMD event generator, version 1.08 [Q. 

2 Experiment 

The NA44 detector is described in [0. The spectrometer is optimized to measure single- 



and two-particle distributions of charged hadrons produced near midrapidity in pA and AA 
collisions. The momentum acceptance is ± 20% of the nominal momentum setting, and only 
particles of a fixed charge sign are detected for a given spectrometer setting. 

The data used for this analysis were taken at the 6 GeV/c setting with the spectrometer 
axis at 44 mrad with respect to the incident beam. The corresponding acceptance for protons 
is 2.35 <y < 2.70 and 0.0 < pr < 0.66 GeV/c {{pr) ~ 230 MeV/c). Scintillator hodoscopes 
are used for tracking, and also provide time-of-flight information with utof < 100 ps. The 
momentum resolution 6p/p is approximately 0.2%. Two threshold Cerenkov detectors in 
conjunction with the time-of-flight (mass-squared) information provide particle identification. 
Particle contamination is less than 1% for the data shown here. An interaction trigger is 
provided by two rectangular scintillator paddles sitting downstream of the target and covering 
approximately 1.3 < rj < 3.5. Offline, a silicon pad multiplicity detector with 27r azimuthal 
coverage over the range 1.5 < rj < 3.3 is used for event characterization. 

Target thicknesses of 10mm and 5mm were used for the p+Pb and S+Pb data, correspond- 
ing to approximately 5.9% and 5.0% interaction lengths respectively. Multiple scattering in 
the targets, combined with the spectrometer momentum resolution, lead to a resolution in k*, 
the particle momentum in the rest frame of the pair, of a{Sk*) ^ 9 MeV/c. The experimental 
trigger required an interaction in the target, at least two tracks in the spectrometer, and a veto 
on pions (and lighter particles) in the spectrometer acceptance. No centrality selection was 
made in the trigger, however the two-track requirement itself biases the data towards high- 
multiplicity events. Offline track reconstruction was followed by track quality cuts, including 
rejection of close-by tracks, and particle identification cuts to select proton pairs. Additional 



cuts to reject multiple interactions in the target (pileup) and require clean events were made. 
The final event samples are approximately 8k pairs for the p + Pb data and 15k pairs for the 
5" + Pb interactions. 

3 Correlation Analysis 

The experimental correlation function is constructed from identified two-proton events. Due 
to limited statistics, only a one- dimensional analysis has been performed. The correlation 
function is calculated as: 

c(n = p^. (1) 

-'■ ^ uncorr x^^ ) 



where k* is the particle momentum in the rest frame of the pair: k* = -^ — , Q = {qo, q} = 
Pi—P2 andpi,p2 are the particle 4-momenta. The correlated, Ncorr, and uncorrelated, Nuncom 
proton pairs are taken from the same and different events respectively. The stability of the 
resulting correlation function under the applied cuts is studied by varying each cut (such as 
the mass-squared selection, and the hodoscope multiplicity) individually. In all cases, the 
correlation function data points are within one standard deviation of their nominal values. 
The stability with respect to different sampling procedures for producing the mixed-event 
distribution is also tested, and introduces no systematic deviations. 

The experimental correlation functions presented here do not contain any correction fac- 
tors. Such corrections are often made in boson correlation measurements to take account of 
the effects of finite detector acceptance and momentum resolution, residual correlations in 
the background distribution |2^, and the Coulomb interaction, thereby recovering the "ideal" 



correlation function in which correlations are due purely to quantum statistics. Rather than 
correct the experimental distributions, theoretical correlation functions based on model pre- 
dictions have been generated, and modified to include each of these effects. The results of 
these simulations are then compared directly to the measured correlation functions. 

The theoretical correlation function is generated by selecting proton pairs from a given 
source model, calculating the associated weight due to quantum statistics and final state 
interactions (strong and Coulomb) 0, |^, and propagating the particles through a realistic 
simulation of the experimental detector, including all instrumental and acceptance effects. A 
full calculation of Coulomb effects is used rather than the Gamow factor approximation. Data 
cuts are then made in exactly the same way as for the experimental data. For this analysis, 
two models have been used as inputs for the simulation: a Gaussian source model and the 
RQMD event generator. 

In the case of the Gaussian model, a Gaussian density distribution p(r, t) = ,^^.2.^3^- 
is assumed for both the spatial and temporal extent of the source [^ . It is assumed that pro- 
tons are emitted by sources moving in the longitudinal direction with the velocities of the 
proton pairs (the Longitudinally Co- moving Source, LCMS). In this way the fast longitudinal 



r2 _ t2 
571 2^ 



motion of proton sources at mid-rapidity is taken into account. Because the rapidity accep- 
tance in this measurement is narrow (2.35 < y < 2.70), the results assuming the LCMS should 
not differ significantly from the results assuming a fixed frame of reference centered at about 
y = 2.5. Since a one-dimensional analysis is performed here, a spherically-symmetric source 
parametrized by equal Vo radius and r time parameters is assumed. The RMS radius in this 
parametrization is given by ^/Srg. 

In order to compare with RQMD predictions, the theoretical correlation function for pairs 
of protons from generated events is calculated 0, and the appropriate weighting for different 
impact parameter collisions is taken into account by matching the multiplicity distributions 
of the data and model. This matching shows that, in contrast to the p + Pb data where only 
the highest multiplicity events are selected by the trigger, the S + Pb data are only weakly 
biased toward central collisions. 

In calculating correlation functions for comparison to data, three factors are particularly 
important: the admixture of indirect protons coming from hyperon (mainly A — >■ n~p) decays, 
the acceptance and resolution of the experiment, and the residual correlations arising in the 
single particle background distributions. 

Weak decays of strange baryons are a significant source of protons and contribute to the 
yields measured in the NA44 spectrometer. The influence of the admixture of indirect protons 
on the shape of correlation function has been studied using data from the RQMD and Venus 
(v5.21) 1^ event generators, combined with a detailed simulation (GEANT) of the detector. 



The two models give similar results: about 22% of protons measured in the spectrometer 
come from the decay of A's in both p + Pb and S + Pb collisions, and cannot be distinguished 
from direct protons. This contribution has only a weak pt dependence, and is greater at 
low Pt- In order to take into account this non-correlated contribution to coincident pairs, 
a fraction of 22% of "decay" protons is included when calculating the correlation function, 
giving about 39% of uncorrelated pairs in the two-particle sample. These pairs are assigned a 
weight of unity in the correlation function calculation. This contribution significantly reduces 
the magnitude of the observed correlation effect, but does not change the general shape of the 
correlation function. 

The acceptance and resolution of the NA44 spectrometer for protons is estimated using 
a detailed simulation of the detector, including systematic effects introduced by the track re- 
construction procedure for close-by pairs of particles. The k* acceptance for proton pairs 
is convoluted with the detector resolution and the rather complicated shape of the two- 
proton correlation function, and gives rise to an asymmetric distortion of the correlation 
peak. This "smearing" is particularly important for small k*, where the expected depletion 
due to Coulomb and statistical effects is masked by the experimental resolution. 

Residual correlations arise in the mixing procedure because the single particle distributions 
forming the reference (background) sample are in fact projections of the measured two-particle 
distribution, which includes the correlations. This effect is taken into account by calculating 
and applying weights which correspond to the mean value of the correlation function for all 
combinations of a given particle with all others [^]. This residual correlation modifies slightly 
the height of the correlation peak, and induces a slope for large values of k*. 



All three effects are included in the simulation procedure. In Fig. 1 they are demonstrated 
separately in order to illustrate their particular features and the relative significance of each 
of them. 



4 Results 

The experimental proton-proton correlation functions for p + Pb and S + Pb collisions are 
presented in Fig. 2. The correlation functions have been normalized such that C{k*) = 1 for 
80 < A;* < 160 MeV/c. The peak observed in the region of small k* values (~ 20 MeV/c) can 
be attributed to the combined effect of the final state interaction and Fermi-Dirac statistics. It 
reflects the space-time properties of the proton emission process. The two correlation functions 
have qualitatively the same shape, but the peak is much more pronounced for the p + Pb data 
than for the 5* + Pb, indicating that protons are emitted from a smaller source in the case of 
p + Pb collisions. Also shown are the correlation functions and associated source sizes from 
the Gaussian model obtained with a minimum-x^ fit over the range < k* < 160 MeV/c. 
The fits yield r^ = r = 1.42t[|:[j^ fm {x^/Ndf = 10.3/10) for p + Pb and r^ = t = 2.651^^ 
fm {x^/Ndp = 9.6/10) for the S + Pb data. (Note that the Tq values have to be multiplied by 
v3 to obtain the corresponding RMS radii.) Fits assuming r = give Tq parameters which 
are equal within errors to those obtained with r = Tq. This is a consequence of the weak time 
dependence of the final state interaction effects for small particle velocities ; for these data 
the mean velocity of the proton pairs in the LCMS frame is ~ 0.24c. 

These results have been obtained assuming that the fraction of measured protons from A 
decay is that given by the simulations described above (22%). While the RQMD and Venus 
models give consistent results for the fraction of protons from hyperon decays in the NA44 
acceptance, there is a factor of order two discrepancy between the available experimental 
data on A production yields in S + A collisions at CERN SPS energies [|6|, |3- The RQMD 
prediction is roughly midway between the two measurements. Using the experimental data as 
an estimate of the uncertainty on the A yield (approx. ± 40%), and hence on the fraction of 
protons from A decays, gives Tq = r = 1.421q;{7 fm for p + Pb and Tq = r = 2.65to;i7 fm for 
S + Pb. The RMS radii of the proton and the sulphur and lead nuclei are 0.88 fm, 3.26 fm 
and 5.50 fm respectively [^. In the frame of this geometrical interpretation, the measured 



source sizes are greater than the size of the projectile but smaller than that of the target. 

In order to analyze the centrality dependence of the correlation effect, both the p + Pb 
and S + Pb data have been divided into three multiplicity bins, based on the charged-particle 
multiplicity measured by the silicon detector. The binning is chosen to give roughly equal 
proton-pair statistics in each sub-sample. The correlation functions corresponding to the three 
multiplicity selections are shown in Fig. 3 ioi p + Pb and S + Pb. A systematic decrease of 
the correlation effect with increasing multiplicity is observed, corresponding to an increase in 
the size of the emitting volume. The Tq values indicate the scale of these changes. 

The comparison of the data with the predictions of RQMD is shown in Fig. 4. In the case 
of p + Pb collisions, the experimental data are well reproduced by the model. For 5* + Pb, 



the model gives a somewhat larger correlation effect than is observed in the data. Fitting the 
S + Pb RQMD prediction with the same Gaussian model (and same fraction of protons from 
A decay) as used for the data yields ro = t = 2.27!!!o;o2 fm, compared to Tq = r = 2.65io;i7 for 
the data. 

5 Discussion 

A simple, one-dimensional description with Gaussian co-moving sources is used to quantify 
the correlation effects. In introducing the co-moving source, it is assumed that particles 
with similar longitudinal velocities are emitted from nearby space-time points, such as occurs 
when fast longitudinal expansion takes place. In this case, the correlation function effectively 
measures only a part of the space-time extent of the emission volume. For the S + Pb system, 
the multiplicity cuts select different impact parameter collisions, and hence different sizes of 
the overlapping region between projectile and target nucleons. This is reflected in the observed 
centrality dependence of the correlation functions. Comparison of the measured source sizes 
with the RMS radii of the projectile and the target is also in qualitative agreement with such 
an interpretation. Note that in the target fragmentation region the sizes measured by two- 



proton correlations reflect rather the size of the target nucleus [^. In the case oi p + Pb 
collisions, the number of projectile participants is always equal to one and the measured size 
can be related to the extent of the hot region along the path of the proton in the lead nucleus. 
The source sizes in this case are much smaller than in the case of the sulphur projectile. 

These data can be compared to correlation results obtained for other particle types near 
midrapidity and in similar pt intervals, which assume the same type of emitting source. For 
both systems studied {p + Pb, S + Pb) , the proton source sizes appear to be smaller than those 
of pions, especially in the case of the sulphur projectile: for p + Pb, rrs(7r"*") = 2.00 ± 0.25 fm 



while for S + Pb, rTs(vr+) = 4.15 ± 0.27 fm {{pr) ^ 150 MeV/c in both cases) p], ig. The 
To values for protons are similar to the corresponding kaon source sizes: for p + Pb, rTs{K~^) 
= 1.22 ± 0.76 fm while for S + Pb, rTs{K+) = 2.55 ± 0.20 fm {{pr) ^ 240 MeV/c in both 
cases) |T^ . (Note that for the three-dimensional boson source parameters, multiplication by a 
factor of vS is also necessary to obtain RMS sizes). The Gaussian source model used in these 
analyses describes sources for which the momenta and positions of the emitted particles are 
independent in the LCMS. Collision dynamics may, however, induce correlations between the 
measured source size and the transverse momentum of the particles. In particular, a decrease 



of the measured size with increasing transverse mass rriT (where rriT = Jpx + "^^ ) ? observed 
in common for pions and kaons, may indicate the collective expansion of an equilibrated 
system formed during the collision jn| |TU|. The smaller source sizes for the two-proton data 



{{mx) ~ 970 MeV) compared to the two-boson data [|T8], |1^, ^ are in qualitative agreement 
with such an itit dependence. 

The contribution of indirect protons coming from hyperon decays is also a parameter in 
the description of the experimental results, and can influence the values of the extracted 
source sizes. The experimental uncertainty on the A yield is described above, and results 



in a systematic error of approximately 0.2 fm on the extracted To parameters. Of the other 
hyperons, Monte Carlo studies show that only the S+ also contributes to the proton yield 
measured in the NA44 spectrometer. RQMD and Venus predict that the number of detected 
protons from E+ decay is at most 30% and 10%, respectively, of that from A decay - less than 
the experimental uncertainty on the A yield itself. Possible uncertainty in the contribution 
of indirect protons cannot significantly alter the conclusion on the size of the proton source 
compared to the pion source however, nor on the centrality dependence of the observed sizes. 
In order to produce the change in the measured correlation functions between the lowest and 
highest multiplicity bins assuming that the difference is due entirely to indirect protons implies 
that more than 50% of detected protons come from hyperon decays. However, an enhanced 
emission of hyperons with respect to that predicted by the RQMD model in the case of S + Pb 
collisions may be responsible for the difference between the model prediction and the results 
of measurements. 

6 Conclusion 

These results demonstrate the first attempt to include baryon-baryon correlations in the anal- 
ysis of the space-time dynamics of particle emission at midrapidity at CERN SPS energies. 
The measured source sizes increase with the mass of the projectile and with the collision mul- 
tiplicity. Within the context of a Gaussian source model, the proton radius parameters are 
smaller than the size of the target nucleus but larger than the sizes of the projectiles. Com- 
paring to other particle types, the proton size parameters are smaller than those for pions, and 
similar to those for kaons, measured in the same collisions. (Note that collision dynamics may 
induce a dependence of the measured source sizes on the momentum of the emitted particles.) 
Good agreement with the results of RQMD (vl.08) simulations is seen for p + Pb collisions. 
For S + Pb collisions, the measured correlation effect is somewhat weaker than that predicted 
by the model simulations. 
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Figure 1: The influence of different experimental factors on the form of the measured corre- 
lation function. The example of a Gaussian source with r^ = t = 2.38 fm is shown, where the 
effect of each contribution has been added cumulatively. 
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Figure 2: Two-proton correlation function for p + P6 (left) and S + P6 (right) collisions. 
The open points and dashed lines show the predicted correlation functions for longitudinally 
co-moving Gaussian sources with the radii indicated in the figure. 
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Figure 3: Two-proton correlation function for p + Pb (left) and S* + Pb (right) collisions 
for three different multiplicity selections. The open points and dashed lines show the form of 
correlation function for longitudinally co-moving Gaussian sources with the radii indicated in 
the figure. 
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Figure 4: Comparison of the experimental correlation functions with predictions of the 
RQMD (vl.08) model for p + Pb (left) and 5* + Pb (right) collisions. The hatched error 
band on the RQMD predictions includes the systematic error due to the uncertainty in the 
spectrometer resolution and the uncertainty on the impact parameter weighting obtained by 
matching the multiplicity distributions of the data and the model. The non-hatched error 
band reflects the experimental uncertainty on the measured A yields [^, 12^. Statistical 
errors are negligible. 
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